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In this paper we give an approximation method for obtaining the 
probability the server is busy and the mean waiting time as seen by 
the nth arriving customer for the GI/G/1 queueing system. Transient 
behavior is the key issue of the method. The approximation consists 
of a pair of recursion formulae whose state variables are the proba- 
bility of delay and the mean waiting time. Any initial state may be 
prescribed for the Oth arriving customer. Programming is very easy, 
and the computation is rapid. The procedure is useful for rush-hour 
analyses and for studying the recovery of a system from temporary 
overload. 

I. INTRODUCTION 

This paper presents an approximation method for obtaining the 
probability the server is busy and the mean waiting time as seen by 
the «th arriving customer for the GI/G/1 queueing system. Thus, 
transient behavior is the key issue of the method. The approximation 
consists of a simultaneous pair of recursion formulae whose state 
variables are the probability of delay and the mean waiting time. It is 
assumed that the Oth arriving customer finds the queue in some 
prescribed state from which the successive states are computed. Nat- 
urally, for the nth arrival when n is large, the computations provide 
approximations to the corresponding eqilibrium quantities when equi- 
librium exists. The procedure, however, is not limited to queues 
possessing an equilibrium state. For methods specially adapted to 
approximating the equilibrium quantities, we refer to a paper by A. A. 
Fredericks, 1 in which the approximation A ,i of that paper essentially 
corresponds to the equilibrium results obtained here; we also refer to 
Fredericks for an application to computer systems. 2 

Transient analysis is of particular importance in studying the recov- 
ery of a system from temporary overload. This can occur after a short 
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downtime during which the buffers become full and the machine must 
now work off the accumulated load. Another situation calling for 
transient analysis occurs in the case of minicomputers receiving hourly 
data from an electronic switching system (ESS). The hourly rush of 
traffic requires a transient analysis to determine the buildup and falloff 
of delays. This corresponds to the general problem of rush-hour 
situations. Knowledge of the transient behavior of a queue also indi- 
cates the number of customers needed in the arrival stream until 
approximate equilibrium conditions are recovered from some tempo- 
rary overload (i.e., the relaxation time may be estimated). 

In Section II we obtain the general recursion formulae and the 
corresponding formulae defining the equilibrium results. In Section III 
we specialize the recursions to the GI/M/1 queue and present the 
exact solutions obtained using the Takacs 3 method as a control. We 
show the relationship between the exact probability of delay and mean 
waiting time, and also conclude that the approximate equilibrium 
quantities are, in fact, exact. We discuss a number of numerical 
examples and compare them with exact results. In Section IV we 
reduce the general recursions to the M/G/l case. A relation is again 
obtained between the probability of delay and the mean waiting time. 
The interesting property is proved that the approximate probability of 
delay and mean waiting time as seen by the nth arrival satisfies the 
same relation as the exact quantities; thus, the approximation method 
is shown to preserve certain properties of the exact solution. We 
discuss a number of numerical examples and compare them with 
exact results. In Section V we present several numerical examples of 
GI/G/1 queues. Exact solutions, however, are derived only for equi- 
librium values. 

In all the numerical examples given, the queues start empty; how- 
ever, this is done only to facilitate obtaining the controls. Any initial 
state may, of course, be used in the recursion formulae. A synopsis of 
the important recursions and definitions of symbols is given in Section 
VI. 

II. RECURSION FORMULAE 

Let A(x) be the distribution of time between arrivals with mean 
arrival rate A, and B(x) the distribution of service time with mean 
service rate ju; the random variable £, which is service time minus 
interarrival time, has the distribution 

K(x)= J B(x + y)dA(y). (1) 



If W n is the waiting time of the nth arrival, then we have the recursion 
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W a = [Wn-l + £n] + , (2) 

in which 

[x] + = x, x>0 

= 0, x < 

and the £„ are iid with the common distribution K(x). Let W n (x) be 
the distribution function of W„; then W n (s) is defined by 



W n (s) = e- sx dW n (x), (3) 

Jo_ 

and hence 

W n (s) = Ee- sW ". (4) 

A recursion relating W n (s) to W n -i(s) will now be developed. From 
eqs. (2) and (4) we have 

W n (s) = Ee- s(W »-' + i" )+ , (5) 

and hence 

W n (s) = E J e- siW -> +x) *dK(x), (6) 

in which the expectation is over the distribution of W n -\. Let 

K(s) = J e- sx dK(x); (7) 

then the functions 

K + (s) = I e~ sx dK(x) (8) 

-0_ 

X_(s) = I e- sx dK(x) (9) 



and 



will be used in the development. From eq. (6), we find that 

W n (s) = Wn^K+is) + I Ee- slw -> +xl+ dK(x). (10) 



The recursion of eq. (10) is exact but, to obtain a simple, explicit 
recursion, an approximation will be used at this point for the distri- 
bution of W n -i- The more accurate the approximation chosen, the 
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more accurate the final recursion will be, but also, presumably, the 
more difficult to use; accordingly, a simple exponential approximation 
will be used. 
Let 

a n = EW n , (ID 

and 

Jn = l-W n (0+); (12) 

then we use the approximation 



_j 2 
W n {x)<*l-J n e *■*, *>0, 



= 0, x < 0. 



From eq. (13) we compute 



Jn-X 



Ee-sWn-^r = ! _ f n ~ lS «,«,_, " x < . 



e/n-l 



<Xn-l 



+ S 



Thus, eqs. (10) and (14) yield the recursion 

Jn-lS 



W n (s) = W n -l(s)£ + (s) + £-(0) - 



e/n-l 



a„-i 



ft- - 



+ s 



e/n-l 



an-i 



(13) 



(14) 



(15) 



A pair of recursion relations will now be obtained for J n , a n . Since 

lim W n (s) = 1 - J nt (16) 



and 



we obtain from eq. (15) 



-W' n (0) - a n , 
lim £ + (s) = 0, 



Jn = K + (0) + Jn-lK- ~ 



e/n-l 



an-i 



and 



a n = 



K + (0) + K-1- 



e/n-l 



«n-l 



On-l - £;<o). 



(17) 
(18) 

(19) 
(20) 



Equations (19) and (20) constitute the approximate recursions sought. 
If the queueing system possesses equilibrium values, then approxi- 
mations designated by J, a, W(s) are obtained from 



J=£ + (0)+e/£-(-^, 



(21) 
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a = 



A\(0) +K [-- 



a - K'+(0), 



and 



W(s) = 



1 - K + (s) 



- + s x 
a 



(22) 



(23) 



The recursions of eqs. (15), (19), and (20), and all subsequent 
recursions derived from them, apply to arbitrary initial conditions and 
to stable or unstable queues. 

III. G//M/1 QUEUEING SYSTEM 

To ascertain the quality of the approximations obtained by eqs. (19) 
and (20), a control is needed. This will be provided by the explicit 
solution formulated by Takacs 3 for the GI/G/1 system. Accordingly, 
let 



G(z, s) = E W n (s)z n , 

n-0 



A(s) = I e 
'o 



x dA{x), 



and 



then 



B(s) 



-f 



e- 8X dB(x); 



&(s) = A(-s)B(s), 



and the factorization 



(24) 



(25) 



(26) 



(27) 



(28) 



1 - zK(s) = r+Cz, s)T-(z, s) 

defines the functions T+(z, s), T-(z, s). The function T + (z, s) is to be 
analytic for R e s > and not to vanish for R e s > 0, while T-(z, s) is to 
be analytic for R e s < and not to vanish for R e s < 0. Further, define 
the projection operator T by 

TEe~ 8v = Ee~ sv+ , (29) 

in which tj is an arbitrarily given random variable. Then the exact 
solution for the generating function, G(x, s), is 



G(z, s) = 



r + (z, s) 



W (s) 



T-(z, s) 



(30) 
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The only exact solutions to be studied in this paper will correspond 
to the queueing system starting empty; hence, W (s) = 1. For this case 
the projection of eq. (30) may be evaluated; we find that 

G < 2 > s > = m — L , m - (3D 

T + (z, s)T-(z, 0) 

To apply eq. (31) to the GI/M/1, we consider 

R(s)=A(-s)^— (32) 

ft + s 

and 

vi n + s — nzA(-s) 



H + s 




Let 5(2) > be defined by 




H-8-ilzA(8) = 0; 


(34) 


then we find that 




r>/ x 8 + s n + s - fizA(-s) 

1 - zK{s) = — - — . 

(i + s 8 + s 


(35) 


Thus, 




T + (2,s)=— — , 
\i + s 


(36) 


r_( 2 , s) =^ s ;^ ( - s) , 


(37) 



and 



r_(z, 0) = ^ (1 5 Z) . (38) 



Hence, from eq. (31), we obtain 

<**•*> — TT-^inr- < 39 > 

H(l — z) 8 + s 

Let j n , a n designate the exact values of 1 — W„(0+), and EW n , 
respectively, and let 

J(z) = i jnZ\ (40) 



a(z) = I a n z n ; (41) 

n-0 

then, from 
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M - 



1 - z 



- G(z, oo), 



a{z) - - — G(z, s)| 



and eq. (39), we get 



and 



/<*) = 



a(z) = 



1 -2 ' 
1-2 



(42) 
(43) 

(44) 
(45) 



Let L n =jn -jn-i,A n - a„ - a„-i and L(2), A(2) be the corresponding 
generating functions; then elimination of 8 in eqs. (44) and (45) yields 



At i : L(2) 
A(2) =-• 



/i 1 - L(2) 



and 



a(z) =-■ 



(46) 



(47) 



M 1 - L(z) 
In particular eq. (47) implies the recursion relating j n , a n , as follows: 

(48) 



1 . n_1 
a n = -j n + X A„-k(jk —jk-i)- 



From eq. (32), we find that 
K + (s) = 



P 



H + s 



A(n), 



and 



Z-(s) =-?—[A(-s) - A(n)l 
H + s 



K + (0)=A(fi), K' + (0) = --A(ii). 



Hence, eqs. (19) and (20) specialized to the GI/M/1 are 



J n =A(n) + 



1 Jn-1 



OLn-l 



(49) 
(50) 

(51) 
(52) 



and 



a n = 



A([i) + 



I Jn-l 

II CLn-1 



Oin-l 



a n -i + -A{p). (53) 
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We will now show that the limiting forms of J n , a n , that is J, a are 
exact; thus J = j, a = a. Equations (52) and (53) show that 

J = A( f i)+ J 






a(^\-A(ji) 



and 



a = 



A([i) + 



1 J 

1 L 

fi a 



A(^)~A(ii) 



Hence, 



a(£)-a w 



Tj 1 -? A( " ) 



1 

/i a 



and, from eqs. (55) and (56), 



1 J 

a = — . 



and 



Hl-J 



J = A[fi(l-J)]. 



(54) 



a + -A([i). (55) 



(56) 



(57) 



(58) 



Since eqs. (57) and (58) are the exact relations defining j, a, the 
statement is proved. 

An M/M/l will be used to start the numerical examples. For the 

M/M/l, one has (p = A//i) 

Jr,=^—l P + ^ \ (59) 



1 + p 



1 + x^ 



and 



Ctn = 



1+P 



P + 



1+jffr-l 



, P 
Oin-\ + ~ 



in which the designation g„ = J n /<x n will henceforth be used. 
Since Jo = 0, ao = 0, the computations are started with 



Ji = K + (0) = 



1 + p' 



and 



(60) 
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Table I— M/M/1, p = 0.2. 



1 


0.1667 


0.1667 


0.1667 


0.1667 


2 


0.1898 


0.2176 


0.1898 


0.2176 


a 


0.1962 


0.2368 


0.1962 


0.2364 


4 


0.1985 


0.2445 


0.1985 


0.2440 


5 


0.1994 


0.2477 


0.1993 


0.2473 



Table II— M/M/1 , p = 0.8. 



1 


0.4444 


04444 


0.4444 


0.4444 


2 


0.5542 


0.7517 


0.5542 


0.7517 


;i 


0.6047 


0.9959 


0.6084 


0.9912 


4 


0.6354 


1.2016 


0.6418 


1.1890 


5 


0.6570 


1.3804 


0.6650 


1.3578 



Table III— M/M/1, p = 2.0. 



0.6667 
0.8148 
0.8719 
0.9012 
0.9191 



0.6667 
1.2593 
1.8233 
2.3727 
2.9131 



0.6667 
0.8148 
0.8807 
0.9172 
0.9362 



a 



0.6667 
1.2593 
1.8189 
2.3603 
2.8922 



1 P 

<*! = -#;(())=- 



(61) 



Mi + p' 

which are exact. The value ju = 1 will be used in all examples. Tables 
I through III below present approximate and exact values for p = 0.2, 
0.8, and 2.0. 

For the next example, the renewal stream is 

A(s) = 



The recursion relations are 

Jn = — 1 + Jn-\ 

b 
and 

i r 
i + 



(1 + s)(l + 2s) 
5 + 2g n -, 



1 



(1 + g„-i)(l + 2g n - i ) 
5 + 2#„-i 



(1 +£„_!)(! + 2#„-,) 



«„->+-. 



Some numerical results are given in Table IV. 

As a further example the case D/M/l is considered. We have 



A(s) = e 



(62) 



(63) 



(64) 



(65) 



The recursion relations are 
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Table IV— G//M/1 , p = 1 /3. 



1 


0.1667 


0.1667 


0.1667 


0.1667 


2 


0.1991 


0.2269 


0.1991 


0.2269 


3 


0.2100 


0.2538 


0.2101 


0.2534 


4 


0.2147 


0.2670 


0.2148 


0.2662 


CO 


0.2192 


0.2808 


0.2192 


0.2808 



Table V— D/M/1 , T = 2, p = 0.5. 


n 


J a j 


a 



1 0.1353 0.1353 0.1353 0.1353 

2 0.1720 0.1903 0.1720 0.1903 

3 0.1867 0.2179 0.1868 0.2175 

4 0.1938 0.2331 0.1940 0.2324 

5 0.1977 0.2419 0.1978 0.2410 
oo 0.2032 0.2550 0.2032 0.2550 



e -T8n-\ _ e -T 

Jn = e~ T + Jn-l — , (66) 

1 - gn-1 



and 



'-(<*♦*£=?) 



a„-i + e~\ (67) 



We find that the exact result fory„ is 

Jn=i^-T^e- kT . (68) 

The equilibrium values for T = 2 are J = j = 0.2032, a = a = 0.2550. 
Table V presents the numerical results. 

V. M/G/1 QUEUEING SYSTEM 

For the M/G/1 queue, we have 

£(s)=-l—£(s); (69) 

X — s 

hence 

i-rita- Az^igw (to, 

A — s 
Define 6(2) > by 

A-fi-Az£(S) = 0; (71) 

then 



A — s S — s 
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(72) 



and hence 



r + (2, s) = 



r_(z, s) = 



X-s- XzB(s) 



8-s 



8-s 
X-s 



and 



r.uo) = x . 



Accordingly, the generating function, G(z, s), is 

X 8-s 

G(z, s) = -- 



8 X-s- XzB(s) ' 

thus, the generating functions j (2), a(z) are 

1 X 



Hz) = 



and 



a(z) = 



1- z 5' 
1 /l 11- zp 



1- z\8 X 1-2 
Elimination of 5 in eqs. (77) and (78) yields 

zp 



a(z) =-- 

A 1 — 2 



1-2 



-m 



which implies that the relation between j n , a n is 

11. 

JU A 

n 1 " . 

a n = t- L Jk. 

Ii X k =i 

This may be compared with eq. (48) for the GI/M/1. 
From eq. (69), we have 

JL(s)=-^-£(A), 

X-s 



K + (s)=t [B(s)-B(X)l 

X — s 



and 



K + (0) = 1 - B(X), 



K' + (0)=\[1-B(X)]--. 

A JU 



(73) 
(74) 

(75) 

(76) 

(77) 
(78) 

(79) 

(80) 
(81) 

(82) 

(83) 
(84) 

(85) 
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Thus, the recursion formulae eqs. (19) and (20) become, for M/G/l, 

5(A) 



J n =l-£(\)+Jn-l 



(86) 



1 + jfr-i 



1 - B(X) + 



B(X) 



1 + x*-\ 



«n-i + - - T [1 - 5(A)]. (87) 



Insofar as approximations imitate characteristics of an original, we 
may better apply the approximations. We will now show that the 
approximations J n , a ni satisfy eq. (80). From eq. (86) we find that 

£(A) J n - 1 + B(\) 



! + -*„-> 



J/i-1 



1 Jn-l\Jn ~ 1 + B(X)] 

Otn-1 = T 



and, from eq. (87), 

dn — OLn-1 = — 



Xj n _ 1 S(X)-J n +l-^(X)' 
Jn-l£(X) -Jn+1- B(X) 



J/i-l 

■a n -i + -- T [l-B(X)l 

JU A 



(88) 



(89) 



(90) 



Thus, substitution of a n -i from eq. (89) into the dexter of eq. (90) yields 

a n = a n -\ H 7 J n , (91) 

jU A 

and 



a n = olq + 2j ^a- 

ju, A *=o 



From eqs. (81) and (92), we have 



1 " 



a n — a„ = -a + Jo — 7 E O* - ^*) 5 

A *=i 

hence, for a stable queue, 

jk - Jk -> 0, A: -> 00. 



(92) 



(93) 



(94) 



Since/* — > p, then also J k — » p, £ — > °°. Unlike GI/M/1, however, a ^ 
a. The equality that occurs in GI/M/1 was to be expected since the 
waiting-time distribution was approximated by an exponential in a 
portion of the integration producing eq. (15), and the exact eqiiilibrium 
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distribution is, in fact, exponential. This does not occur in M/G/l, 

however, since the waiting time is not exponential. The value for a is 

simply 

\B"(0) 

a = -- , (95) 

2 1 - p 

while solution of the equilibrium form of eqs. (86) and (87) shows that 



a = - 



1 p - 1 + B(X) 



jMl-p)[l-B(A)]' 

thus, a a a for small A. 

For the first numerical example, consider 



(96) 



A(s)=-?—, B(s)=- — -2; (97) 



p + s (I 

1 + 2 S 



then 



Jn = P^ [4 + p + J^ — ^— , (98) 

(2 + p) 2 \ ^ P + gn-lj 

a " = (2^( 4 + p+ 7^) an - 1 + 1 "A U + P) 



4 + p+— «n-l+l~- 1-77T— T2 • 09) 



We have 



p(4 + p) _ 1 

Jl= 72T^' ai - 1_ P 



4 
1-- v! 



(100) 



(2 + P Y 

for the queue starting empty. Table VI presents some numerical values 
for p = 0.5. 

As another example let us consider the following M/D/l: 

^( S )=TT-7T' B{s)=e- S . (101) 

1 + 2s 

The recursion formulae are 

^-^-'"-'rSb (102) 

and 



a, - 1 - e~ l/2 +-^r <*.-, + 2e~ l/2 - 1. (103) 

y 1 + Zgn-l/ 

The values are summarized in Table VII below. 

V. Gf/G/1 QUEUEING SYSTEM 

For the next group of examples, the control will provide only 
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equilibrium values. The following uses an interrupted Poisson arrival 
stream 4 defined by 

3.004s + 0.913216 
s 2 + 4.308s + 0.913216 



Ms) = 2 .' ,™» rr^S77 • (104) 



The service time distribution is given by 



Thus, 



A 3 + 5s 

&(s) - n . o wo . « v - 106 

(1 + 2s) (3 + 2s) 



- (0.913216 - 3.004s) (3 + 5s) Q 

{S) = (0.223586 - s)(4.084414 - s)(l + 2s)(3 + 2s) ' 

~ , v 0.0516472 0.672764 ,,„_ 

^ W = 0.223586 -s + 4.084414 -s ' (1 ° 7) 

~ , v 0.182021 1.266801 M-M1 , 

^ s, =TT2T + TTir- <108) 



and 



£ + (0) = 0.604288, £;(0) = -0.645553. (109) 

The recursion formulae are 



r , 0.0516472 0.672764 . / , v 

J " < mmHH + *'" ■ ' 0.223586 + g „-, + 4.0844 1 4 +g „ '' " ,0) 



and 



Jn-l)' 

L^. nnn 0.0516472 0.672764 \ 

a n - 0.604288 + + 

\ 0.223586 + gn-i 4.084414 + g n -i) 

• a„-i + 0.645553. 

Table VI— M/G/1 , p = 0.5. 



n 


J 


a 


J 


a 


1 


0.3600 


0.2800 


0.3600 


0.2800 


2 


0.4245 


0.4310 


0.4266 


0.4269 


3 


0.4515 


0.5280 


0.4547 


0.5174 


4 


0.4666 


0.5948 


0.4698 


0.5777 


5 


0.4762 


0.6423 


0.4789 


0.6198 


00 


0.5000 


0.7778 


0.5000 


0.7500 




Table VII 


I— M/D/1 , 


p = 0.5. 




n 


J 


a 


J 


a 


1 


0.3935 


0.2131 


0.3935 


0.2130 


2 


0.4443 


0.3244 


0.4482 


0.3166 


3 


0.4655 


0.3933 


0.4701 


0.3764 


4 


0.4773 


0.4387 


0.4812 


0.4140 


5 


0.4846 


0.4694 


0.4876 


0.4388 


00 


0.5000 


0.5415 


0.5000 


0.5000 



(111) 
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Table VIII— G//G/1 , interrupted 
Poisson, p = 0.7. 



n 


J 


a 


1 


0.6043 


0.6456 


2 


0.7123 


1.1509 


3 


0.7498 


1.5762 


4 


0.7703 


1.9470 


5 


0.7842 


2.2769 


6 


0.7947 


2.5743 


7 


0.8031 


2.8452 


00 


0.8824 


8.0160 



Table IX — Approximate and 
exact waiting-time distributions. 



X 


w(x) = 


w(x) = 





0.1176 


0.1183 


1 


0.2190 


0.2212 


2 


0.3014 


0.3026 


3 


0.3728 


0.3717 


4 


0.4363 


0.4323 


5 


0.4933 


0.4863 


6 


0.5448 


0.5350 


7 


0.5911 


0.5789 


8 


0.6330 


0.6187 


9 


0.6706 


0.6547 


10 


0.7045 


0.6872 



Table VIII presents the numerical results. The exact equilibrium 
values are/' - 0.8817, o = 8.547. 

We will use this example to illustrate eq. (23) for the approximate 
distribution function. Using the approximate equilibrium values J, a, 
as in Table VIII, we find that 

w(x) « 1 - 0.894602e- on008 * + 0.0439472e "° 306135 * (n 2 ) 

- 0.0317493e- 0969454 *. 
The exact waiting-time distribution is 

w(x) - 1 - 0.8420579e- 009903631 - 0.03963069e-° 8959713 *. (113) 

The numerical comparison is given in Table IX. 
As another GI/G/1 example, we consider 

* s > = (i + S )a + 2 S) - (U4> 

and 

£(«) = tee- 78 + tee" 1 - 8 *. (115) 

Here the service consists of a step function with two values. The 
recursive equations are 
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t rt , r ™o. r / 1.226734 0.384559\ , v 
J n = 0.157825 + Jn-i (- , 116 

\l + 2g n -l 1+gn-l) 

and 

/ 1.226734 0.384559\ 
a n = 0.157825 + a„_i + 0.068909. (117) 

\ l + 2g„-i l+g n -lj 

We show the numerical values in Table X. The exact equilibrium 
values are/ = 0.1798, a = 0.09266. 

Table X— G//G/1 , two-step 

service time distribution, p = 

1/3. 



n 


J 


a 


1 


0.1578 


0.06891 


2 


0.1741 


0.08688 


3 


0.1782 


0.09278 


4 


0.1795 


0.09485 


5 


0.1800 


0.09558 


00 


0.1802 


0.0960 



For the last example, a D/G/l will be considered. Let 



A(s) = e~ 8T , 8(8) = 1 —- J ; (118) 



1 
1 + 2 S 



then 



e » T . e ~ 2T 2Te~ 2T 

K{s)= i — rv ^ +(s)= t — rv + — r~' (119) 



and 



&-(s) = &(s) - £ + (s). (120) 

The recursion equations are 

Jn = e- 2T (l + 2T) + Jn-lZ-(-gn-l), (121) 

and 

a n = [e- 2T (l + 2T) + &-(-g n -i)]a n -i + e~ 2T (l + T). (122) 

Tables XI through XIII show values for T = 2, 10/9, and 0.5, respec- 
tively. 

To obtain exact equilibrium values for the stable queues, define the 
roots 8i(z), 8 2 (z) by 

1 - m(z) = v£e- 1/2 ™' (r) , (123) 
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Table XI— D/G/1, t = 2, p = 
0.5. 



0.0916 
0.1085 
0.1136 
0.1154 
0.1162 
0.1167 



0.05495 
0.07016 
0.07561 
0.07775 
0.07863 
0.07927 



Table XII— D/G/1 ,7=1 0/9, p 
= 0.9. 



n 


J 


a 


1 


0.3492 


0.2288 


2 


0.4614 


0.3822 


3 


0.5185 


0.5025 


4 


0.5545 


0.6032 


5 


0.5801 


0.6906 


CO 


0.7590 


2.0631 


Table XIII- 


-D/G/1 

2. 


T= 0.5, p = 


n 


J 


a 


1 


0.7358 


0.5518 


2 


0.8875 


1.0716 


3 


0.9362 


1.5822 


4 


0.9566 


2.0893 


5 


0.9672 


2.5945 


10 


0.9850 


5.1101 


20 


0.9928 


10.1252 



then 



and 



fc&U) - 1 = A-W*! 



m = 



1-2 



1 - - S 1 (z)8 2 (z) 
4 



a(z) = 



1 



Siur 1 + 5 2 (2) _1 - 1 



(124) 



(125) 



(126) 



1 - z 

Thus, the equilibrium values are 

>=1- , / 4 5 1 (l)5 2 (l), a = Sid) -1 + 52(1) _1 - 1; (127) 

and, hence, for T = 2, we have j = 0.1164, a = 0.07842, and for T = 
10/9,y' = 0.7587, a = 1.9895. 
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VI. SYNOPSIS 

We list the recursions and definitions of symbols here for ready 
reference. 



6.1 Definitions 



A(x): interarrival time distribution. 
B(x): service time distribution. 
W„(x): waiting-time distribution of nth arrival. 



-I 



K(x)= B(x+y)dA(y). 
Jo- 

X: mean arrival rate, 
/i: mean service rate. 

P = A//*. 



&(s) 


■f. 


e-"dK(x). 


£ + (s) 


■r 


e- sx dK(x). 


K-(s) - 


■c 


e-*dK(x). 


A(S): 


■/: 


e-"dA(x). 


B(s)- 


-r 


e'^dBix). 


W n (s) - 


■ r 


e 8X dW n (x). 



j n : probability nth arrival sees server is busy. 

J„: approximate evaluation ofy'„. 

a n : mean waiting time of nth arrival. 

a n : approximate evaluation of a n . 

gn = Jn/ctn. 

All symbols without the subscript n designate equilibrium values. 
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6.2 Recursion formulae 
6.2.1 General — Transient 

J n = K + (0) + Jn-lK(-gn-l), 

a n = [K + (0) + K-i-gn-^Ctn-x ~ K' + (0), 

and 

W n (s) = W n -1(S)K + (S) + K-(0) -^— Jn-lK-(-gn-l). 

S -T gn-1 



6.2.2 General — Equilibrium 

j=& + (o) + J&-(-g), 

a = [K + (0) + K-(-g)]a - £' + (0), 
and 



w(s) = 



1 - K + (s) 



K-(0) —jR-(-g) 



s + g 



6.2.3 G//M/1 

J n = A(fi) + 



e/n-l 



■[A(gn-l)-Mli)], 



1 - - gn-1 



A(u) + 



-[il(^-i)-iC(/i)] 



1 --gn-1 



a„_i + -A(n), 



J = i$[>(l- J)], 



and 









1 


J 






a 


Ml 


-J 


6.2.4 M/G/1 










</« = !- 


B(A) + J„- 


l 


B(X) 





1 + x*"-' 



a n = 



1-B(X) + 



B(\) 



l + x g -> 



OCn-l + - - T [1 " £(A)], 



J = j = P, 



a = - 



1 p - 1 + £(A) 



M (1 - p)[l - B(A)] 
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and 
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